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Abstract 

The higher dimensional cosmology provides a natural setting to treat, at a clas- 
^) i sical level, the cosmological effects of vacuum energy. Here we discuss two situations 

where starting with an ordinary matter field without any equation of state we end 
up with a Chaplygin type of gas apparently as a consequence of extra dimensions. 
In the second case we study the quintessential phenomena in higher dimensional 
spacetime with the help of a Chaplygin type of matter field. The first case suffers 
from the disqualification that no dimensional reduction occurs, which is, however, 
O , , rectified in the second case. Both the models show the sought after feature of oc- 

I I currence of flip in the rate of expansion. It is observed that with the increase of 

^j^' dimensions the occurrence of flip is delayed for both the models, more in line with 

current observational demands. Interestingly we see that depending on some initial 
I conditions our model admits QCDM, ACDM and also Phantom like evolution within 

J> ' a unified framework. Our solutions are general in nature in the sense that when the 

■ extra dimensions are switched off the known 4D model is recovered. Correspondence 

i to a recent work of Quo etal for a quiessence like model is also found. 

p 

\Q . KEYWORDS : cosmology; higher dimensions; accelerating universe 
O : FAGS : 04.20, 04.50 +h 



X' 1. Introduction 

Distance measurement of type la supernovae as well as cosmic microwave background 
anisotropy measurements in the last decade points to a late accelerated expansion of 
the universe. Gravitational interaction being always attractive in nature this finding, 
corroborated by a number of cosmic probes has been puzzling astronomers for a long 
time and has so far evaded any plausible but consistent explanation based on sound 
physical principles. Briefly stated - proposals include among others flavor oscillations 
of axions [Ij, modification of Einstein-Hilbert action via the additional curvature terms 
in the Lagrangian [2], introduction of an evolving cosmological constant in the field 
equations, the role of inhomogeneity to cause acceleration [3j, Brans-Dicke scalar field 
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and finally the presence of a quintessential type of scalar field giving rise to a large negative 
pressure (dark energy) |4j. 

In this context the authors of the present article have been, of late, struggling with the 
idea of explaining the late acceleration as a higher dimensional (HD) effect [5],|6]. In the 
framework of higher dimensional cosmology we have been able to show, though in a rather 
naive way, that the acceleration can be explained as a consequence of the presence of the 
extra spatial dimensions and this effect has been coined as ^dimension driven^ accelerating 
model. In fact here the effective Friedmann equations contain additional terms coming 
from the extra dimensions which may be viewed as a 'fluid' responsible for the current 
acceleration. So here we attempt to incorporate the phenomenon of acceleration within 
the framework of HD spacetime itself without invoking a mysterious scalar field with large 
negative pressure by hand. On the contrary the origin of the extra fluid responsible for 
the acceleration is geometric in origin having strong physical foundation and more in line 
with the spirit of general relativity as proposed by Einstein [7] and later developed by 
Wesson and his collaborators [8J. Moreover in a recent communication [9] it is argued 
that quantum fluctuation in 4D spacetime do not give rise to dark energy but rather 
a possible source of the dark energy is the fluctuations in the quantum fields including 
quantum gravity inhabiting extra compactified dimensions. So we recently witness a spurt 
of activities relating to higher dimensional spacetime in its attempts to unify gravity with 
other forces of nature, development of varied brane models, induced matter proposal and 
very recently the dimension driven quintessential models. The present investigation is 
primarily motivated by two considerations. While we have plenty of multidimensional 
cosmological models in literature [10] and also some sporadic works of Chaplygin inspired 
brane models [HI [12] we are not much aware of models of similar kind in higher dimensions 
either driven by extra dimensions themselves or by Chaplygin type of gas. 

The present work essentially comprises two parts. We have taken a (d+4) dimensional 
homogeneous spacetime with two scale factors and a perfect fluid as a source field. As- 
suming an ansatz in the form of the deceleration parameter we find an expression of the 
scale factor for isotropic expansion. In the process we get a class of solutions for the mat- 
ter field with interesting physical properties closely mimicing the well known expression 
of a generalised Chaplygin gas suitable for accelerating model. One may look upon the 
resulting matter field as a two-component fluid consisting a cosmological constant and a 
perfect fluid obeying a higher dimensional equation of state. One finds that in the early 
phase the pure matter field predominates resulting in a decelerated expansion followed by 
the predominance of the cosmological constant with the expansion transiting to an accel- 
erating phase, while in the intermediate stages our cosmology interpolates between the 
different phases of the universe. For sake of completeness we also construct an equivalent 
scalar field with its potential energy term to simulate the dynamics of the cosmology with 
the actual matter field. 

In the second part of the program the approach is reversed. Here we start with the 
more conventional approach of taking a Chaplygin type of matter field to start with. 
Our first model suffers from the disqualification that both 4D and HD scale factors are 
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expanding so the final reduction to the current manifestly 4D universe does not work. This 
shortcoming is remedied in the second part. While literature abounds with quintessential 
models with Chaplygin type of gas in 4D we are not aware of results of similar kind 
in HD spacetime. However we have not been able to find solutions in a closed form 
with the system of equations finally reducing to a hypergeometric series. In any case 
certain inferences can always be drawn in the extreme cases and our analysis shows that 
an initially decelerating model transits to an accelerating one as in 4D. An interesting 
result in this section is the fact that the effective equation of state at the late stage of 
evolution contains some additional term coming from extra dimensions. Fixing some 
initial conditions the cosmology evolves as QCDM, ACDM or Phantom type. Though 
not exactly similar this points to the ^k-essence^ type of models which lead to cosmic 
acceleration today for a wide range of initial conditions without fine-tuning and without 
invoking an anthropic argument. 



2. Accelerating Universe - I 

We consider the line element of (d+4)-dimensional spacetime 



ds' = dt' - R' — — + r'de' + r'sin'ed(t)'j - A^bdy^dy" (1) 

where ( a, 6 = 4, 3 + c?) are the extra dimensional spatial coordinates and the 3D and 
extra dimensional scale factors R and A depend on time only and K is the 3D curvature 
and the compact manifold is described by the metric 'jab- For our manifold x x S'^ 
the symmetry group of the spatial section is 0(4) x 0{d+ 1). The stress tensor whose 
form will be dictated by Einstein's equations must have the same invariance leading to 
the energy momentum tensor as [T3] 



Too = p , Tij = -p{t)gij , Tab = -Vd{t)9ah (2) 

where the rest of the components vanish. Here p is the isotropic 3-pressure and pd-, that 
in the extra dimensions. 

The independent field equations for our metric (1) are 



R2 1 A2 D4 



R m A 1.,. ^,RA 
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Here we have three independent equations with five unknowns {p,p,pd, A, R) and hence 
we are at hberty to conveniently fix up two initial conditions. Firstly we take a perfect 
fluid with isotropic pressure in all dimensions such that p = Pd and secondly we assume 
a particular form of deceleration parameter as 

where a, b and w are arbitrary constants. The equation (6) yields 



R{t) = Rosmh^cut (7) 

i'lT' — ^) such that we get 

1 — n cosh^ut 

1 = I^T— (8) 

n cosh uji 

showing that the exponent n determines the evolution of q. While for n > 1, it is only 
accelerating but for n < 1 we are able to achieve the desirable feature of flip, although 
it is not obvious from our analysis at what value of redshift this flip occurs. The second 
assumption p — pa gives 

R R^ A ^^A^ , ^ 

Relevant to point out that R = A identically satisfies this equation and in this case the 
matter density and isotropic pressure come out to be 



p = l(d + 2)(d + 3)|^ (10) 

p 1 p2 

P = P. = -(c^ + 2)- - -(d + + 2)- (11) 

Now we are in a position to calculate the expressions for matter density and pressure and 
find 

p = A + ^ (12) 

XL" 

n(ci + 3)-2 B 
n(o( + 3) Rn 

2^ 

where, K = \{d + 2){d + ?>)v?uj'^ and B = AR^ . 
We also get the equation of state as 

2A 
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where, (7 — 1^ 



iSyTT' «o ^ = and w = 7 J such that we finally get 



R{t) = Ro 



sinh^ 



{d + 3)A 



It 



7(d+3) 



2(d + 2) 

The last equation makes interesting reading. As t ~ , sinht ~ t and we get 

R = t^K^ 

such that for 4D case((i = 0)[Ti] 

R = t^ 



(15) 



(16) 



(17) 



which yields the well known solutions as t = , t = ti for radiation and dust respectively 
for 4D. Replacing n in the matter field equations we finally get 



p = A + 



B 



and p = = -A - (1 - 7) 



B 

i?7(rf+3) 



;i8) 



From our analysis it follows that 
(i) for dust, 7 = 1 and n = 



P = A + -^(^ and p=-A (19) 



[ii) for radiation, 7 = and n = 



1 \ B 
p = A H ---Y and n = -A + — (20) 



'in) for stiff fluid, 7 = 2 and n = 



B B 

p = ^ + ^m^) p = -K + ^^^ (21) 



One can calculate q for different values of 7. Taking 
(?) 7 = 2 (stiff fluid) 

id 

a = — 

1 + sinh oot 



(d + 2) - sinh^wt , , 



'^0 7 = jfl ( radiation) 



_ {d + 2) -2 sinh^cjt 
^ ~ 2 + 2 sinh^wt ^ ' 
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{m) 7 = 1 (dust) 

_ (d + 1) - 2 smh^cot 
^ ~ 2 + 2 sinh^wt 
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Figure 1: The variation of q and t for different values of d in dust dominated era. Flip is 
delayed as dimensions increases. 



At this stage one can check the influence of extra dimension on the occurrence of flip. The 
figure-1 shows that with dimensions the instant of flip is delayed. The fact that the time 
of flip is delayed with dimensions is in line with observational demands. This is because 
there is plenty of observational evidence for a decelerated universe in the recent past [15] . 
However the dominance of vacuum energy over ordinary matter at some z^q > does not 
guarantee the present acceleration of the universe. For this, the vacuum energy has to 
dominate long enough as to overcome the gravitational attraction produced by ordinary 
matter. Current dynamical mass measurements suggest that the matter content of the 
universe adds up to 30 percent of the critical density. It is generally believed that the 
transition from deceleration to acceleration occurs at a redshift 0.28 < zt < 0.68, which 
confirms the idea that the acceleration is a recent phenomenon. Hence the significance of 
the above plot is particularly encouraging in the context of higher dimensional cosmology. 

Now the expression of speed of sound will be 

C'i-'i = --^ (25) 

dp 7 

which gives (i) for dust case: = 0; (u) for radiation case : = and {Hi) for 
stiff fluid : = 1. Thus we see that the velocity of sound vanishes in the dust model as 
expected and it is maximum in 4D and decreases with dimension in radiation dominated 
era. To avoid imaginary value of the speed of the sound 7 < 2. Again Cs should never 
exceed the speed of light for the range of 7, 1 < 7 < 2. 

As is well known the whole dynamics can be simulated in a field theoretic approach with 
the help of an equivalent scalar field [TB] . For a scalar field (p having self-interacting po- 
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tential U{(j)), the Lagranian of the scalar field will be £ = |(9^(/)9^0 — f/(0), the equivalent 
energy density and pressure for the scalar field will be 

= + = A + (26) 
and the corresponding 'pressure' as 

= - U(^) = - + 7^Jt3) (27) 

such that 

= (28) 
which, in turn, gives via equation (3) for (d+4) spacetime 



V 2 R-JAW'M3) + B ^ ' 

(here a prime and a dot overhead denote differentiation w.r.t. R and t respectively.) 



or 



and 



tanh-^ (coshwt) (31) 

7(a + 3) 



We also see from figure-2 that for dust and radiation cases U {(f)) decreases with (p unlike 
the stiff fluid case where it remains constant. 

It may not be out of place to call attention to a quintessential model driven by a tachyonic 
scalar field [IT] with a potential in 4D space time 



V{T) = ,^ , Wl-7cos2 )T (33) 



,in2 [h/f^ T 

(T is a tachyonic scalar field) giving the cosmological evolution as 



_2_ 

37 



i?(t) = i?o I sinh^^^ I (34) 
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- 7 = 2 (stiff Fliud) 


- 


y = 5/4 (Radiation) 
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Figure 2: The variation of U{(f)) and (p is shown in this figure for different values of 7 in 5D 
case. 



It behaves like a two fluid model where one of the fluids is a cosmological constant while 
the other obeys a state equation p = (7 — l)p, (0 < 7 < 2). Similarity of this evolution 
with our model is more than apparent except for some numerical factors coming out 
because we are here dealing with a higher dimensional spacetime. But the main result 
may be re-emphasized that we get this evolution without forcing ourselves to invoke any 
extraneous tachyonic type of scalar field. 

3. Accelerating Universe - II 

As commented in the introduction the first model suffers from the shortcoming that no 
dimensional reduction occurs. We here adopt the conventional approach of taking a 
predetermined form of energy momentum tensor and then try to solve the equations of 
motion with the gauge condition 

A{t) = Rit)-"" (35) 



where m is any positive number so that dimensional reduction is ensured a priori. For 
the matter field we here assume an equation of state of a generalised Chaplygin type of 
gas in 3D space only 

p=(7-l)p-i?p-° (36) 
where < a < 1 along with a higher dimensional isotropic pressure pd existing in the 
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extra space such the the field equations reduce to 



R 1 

-p= {2 - dm) - + - [m^d{d + 1) + 2{1 - dm)]— . 

R 2 R 

R 1 it^ 

-Pd = (3 - dm + m)- + -[m(d - l)(dm - 4) + 6] — 



(37) 
(38) 
(39) 



where k = dm^[d — 1) + 6(1 — dm). 



For positive energy density, k must be greater than zero which gives m < ^'^ d(d-^^ ~^ 



m > 



d(d-l) ■ 



Again Bianchi identity will be in our case 



R A 
P + 3^{p + p) + dj{p + Pd) = 



(40) 



Now using equations (36), (37) & (40) we get 



P + 



k 



R 



(2 - dm) R 
Solving equation (41) we get 



2dm(m + 1) , „ _ 
7 + ^ )P-Bp- 



(41) 



P 



Bk 



+ 



Ti/T ■ (l+a)M 
R (2-d"i) 



1 

1+a 



(42) 



where 

M = 'jk + 2dm{m + 1) (43) 
and c is the integration constant. From physical considerations we can finally put the 
restriction on m as m < d{d-iy ~^ ior d ^ 1 and m < 1 for c? = 1. 

The relation (43) needs some interpretation. The last term is a typical higher dimensional 
effect, absent in 4D {d — 0). Now we see that several possibihties for different values of 
m present themselves. 

{i) m = —1: Here A = R and we have to sacrifice the desirable feature of dimensional 
reduction to get an isotropic expansion in all dimensions. However it not exactly the 
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cosmology of the section I because we have a Chaplygin type of gas to start with in this 
case. Moreover the isotropy of the metric dictates that p = Pd and so we end up with 
an isotropic pressure in all dimensions. The solutions closely resemble the earlier work of 
Debnath et al [I8] in 4D. 

[ii] m = 0: Here we get flat extra space although the total number of dimensions 
continues to be (c? + 4). But the cosmology is exactly similar to the 4D case referred to 
earlier [18]. In fact this similarity is a direct consequence of a little known theorem of 
Campbell that any analytic N-dimensional Riemmanian manifold can be locally embedded 
in a higher dimensional Ricci-fiat manifold [19]. 

{in) d = 0: Here we simply recover the 4D metric and all the known solutions of 4D 
follow. 

Now with the help of equations (37) & (42), we get 

^ _ 2 

We have not been able, so far, to find a solution of equation (44) in a closed form. Rather 
a Hypergeometric series solution results given by 

RW^ = + c (45) 



where s = 2{i+a) 2-^1 is the hypergeometric function. 

Even then, fixing the values of different parameters one can get the temporal behaviour 
of the scale factors as given in the adjoining figure-3. A cursory look at the figure shows 
that at a certain stage of evolution the cosmology starts inflating. Another desirable 
feature is the fact that the extra dimensions compactify at very early stage of evolution 
in conformity with both theoretical and observational requirements. 

Now for small value of scale factor R{t) {A{t) should be large in this situation), 

1 _ M 

p ^ ci+"i? (^-d™) (46) 

which is very large and corresponds to the universe dominated by an equation of state, 
p = (7 — l)p as is evident from equation (36). At the late stage of evolution (when R{t) 



Bk 



+ 



{l + a)M 
R (2-dm) 



(44) 



2(2 - dm) 2F1 



s, s, 1 + s. 



i?fcR2''(2-dm) 

7m 
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12 3 4 

Figure 3: The variation of R{t) and A{t) vs t is shown in this figure. Here we see that R (t) 
expands indefinitely and dimensional reduction is possible for extra dimensions. We have taken 
a = 0.5, m = 0.5, d = 1 and 7 = 1. This graph is drawn with the help of 'Mathematica' using 
equations (35) & (44)- 



is large and A(t) small) the expressions for p and p are given below. 

'M 



P 



Bk \ 
~M 



and p 



Bk^ 



M 

T 



(7-1) 



1 + 



k 



-(7-1) 



2dm{m + 1) 



p = Wp (47) 



where 



1 + 



2dm{m + 1) 
k 



(4J 



The last equation is interesting and introduces some new physics in our analysis in the 
context of multidimensional cosmology. Here W 7^ —1. Obviously this is due to the 
presence of extra dimensions in the above equation. As mentioned in the introduction we 
have been working on the idea of ^dimension driven^ acceleration for the past few years 
and the expression (47) points to the fact - how extra dimensions aid the inflationary 
process. 

In 4D case (rf = 0) W = — 1 and a ACDM model results. Otherwise the magnitude of 
W is parameter dependent. When m = 0, i.e. A{t) is a constant we again get back the 
4D case. This is, however, a consequence of the Campbell's theorem [T^ as mentioned 
earlier. 



11 



When m>0,W< —1; So a phantom hke cosmology resuhs with the occurrence of 
'big rip' etc. But the cosmology becomes physically interesting when — 1 < m < such 
that > W > — 1 and we get a quiessence type of model [20]. In this connection one 
should note that predictions from current measurements illustrate a degeneracy in which a 
flat QCDM model in indistinguishable from closed ACDM model, although the predicted 
value of W based on CMB measurement is about W = — |. However, the error bars are 
too broad to reach firm conclusions. 

It will be not out of space to call attention to a recent work by Guo and Zhang et al[2T] 
where a very generalised form of Chaplygin relation is invoked where the constant B is 
assumed to depend on the scale factor i.e., B = B{a). Taking B = B^a'^, for example, it 
is shown that for a very large value of the scale facor the model interpolates between a dust 
dominated phase and a quiessence- dominated phase{i.e. , dark energy with a constant 
equation of state) [22] given hj W = —1 + n/6. So our model closely mimics the above 
work where the extra dimensions takes the role of a variable B{a). 
Now, if we calculate (3+d)-dimensional volume 

V = R^A^ = i?(3-'^™) (49) 

Since m < d(d-r) '~^ ^'^^ d 1 and m < 1 for d = 1, .". the value of (3 — dm) should be 
positive. 

Now for accelerating universe, R> 0, implying 

c M(M + 2dm - 4) 
^ ^ 2B k{2-dm) ^^"^ 

The above expression shows that the universe will be decelerating for small values of scale 
factor while for large values we get accelerating universe and the flip occurs at t = tc such 
that 



c M{M + 2dm-A) 
2B k{2 - dm) 



2 — dm 



M{l + a) 

(51) 



While looking at the plot (figure-4) one must admit that it is rather artificial to take the 
number of spatial dimensions as a continuous variable. However given the fact that in 
Friedmannian cosmology R{t) is always a monotonically increasing function of time, R{t) 
may be taken as a measure of the flipping instant and incidentally it does not initially 
show much variation with d. The fact that it shoots up at some higher value of d need 
not be taken too seriously because it may be due to the arbitrarily chosen initial value of 
a, m etc. In any case the flip is delayed with dimensions and from observational point it 
is encouraging because all the evidences from different cosmic probes point to a very late 
acceleration. 

The speed of sound = |^ in the Chaplygin gas 
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^s' = (7-l)(l + «)-y (52) 



Using equations (38), (39) & (52) we get, 



V 



2 



+ J {m^d{d + 1) + 2(1 -dm)} + ^(2 - dm) (53) 



where q is the deceleration parameter defined as — For accelerating model the value 
of q should be negative. So for real value of Vs, + f {m'^d{d + 1) + 2(1 — dm)} > 

^(2 - dm) which gives + {m'd{d + 1) + 2(1 - dm)} >\ q |. This 

relation gives an upperbound for the q, which may have observational consequences. Since 
< a < 1, it can not exceed the speed of light. But for larger values of a the speed of 
sound exceeds the speed of light. It also depends upon the value of m. In cosmology a 
speed of sound exceeding the speed of light does not contradict the causality [23]. 
Let a = 0.5, m = i, and d = 6, therefore k should be | and M equal to || in this case. 
If we calculate ff, we get two conditions (i) when | g |< 0.923 then Vg < c, (ii) otherwise 

Vs > C. 

Following the analysis given in the previous section we get the analogous energy density 
Pff, and pressure for the scalar field as 



= + f/(0) 



Bk 
~M 



1 

l + a 



+ 



(1^ 



-dm) 



(54) 
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1 • 



</)^- W) = (7-l)p-i?p- 



(7-l) 



Bk 



+ 



/\^ ' (1+q)M 
^^-1- Jl (2-dm) 



1 

1 + a 



B 



Bk 



+ 



^^J- Jl (2-dm) 



(55) 
(56) 



such that we get 



Bk 
~M 



1 k 



(l + a)M 
(2-dm) 



Bk 
~M 



l + Qf 



(l + a)M 
(2-dm) 



5A; 



+ 



C7 



(l + a)M 
(2-dm) 



C7 



(l + a)M 
(2-dm) 



(57) 



(58) 



and 



[/(</)) = -(2 -7) 



5A; 



+ 



M (1+q)M 
Jl (2-dm) 



1 

" 1 + Qf 



B 



Bk 



(l + a)M 
(2-dm) 



(59) 



Integrating equation (58) 



± 




B (7^ - 1) + 



ay 



(l + a)M 
R (2-dm) 



Bk 

M 



(l + a)M 
(2 — dm) 



dR 



(60) 



where 0o is an arbitrary constant of integration. It is very difficult to integrate the 
equation (60) in a closed form. However if at this stage we take the m = —1 i.e., ^ = 1 
then we are able to integrate and get the following solution. 



and 



2 2 + d , if fey {i+a)k-, 

, sinh"^ <^ X / —R 

A;7 1 + a I V ^ 




2-7/5 



+ 



2 V7 

B f B 

7 



coshi+a 



cosh 1 




(61) 



(62) 



If we put c? = 0, i.e., in 4-Dimensional spacetime we get solutions of the above equations 
(61) &; (62) which are identical with 4D case [18]. 
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Figure 5: The variation of (p and R is shown in this figure. The scalar field decays with time 
more sharply in lower dimensions. (Taking a = 0.6, k = 6 and j = 1). 
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Figure 6: The variation of U{(j)) and (j) is shown. It suggests critical 7 dependance. With 
dimensions U{(j)) becomes flatter, (a = 0.6, k = 6). 



4. Discussion 

We have here studied two higher dimensional cosmological models amenable to late ac- 
celeration - one with a particular form of deceleration parameter and for the other with a 
generalised chaplygin type of gas. The important findings may be briefly summarised as: 

{i) Unlike the first case we have dimensional reduction in the second. 

{ii) The striking thing about the first model is the fact that we do not have to assume 
an extraneous, unphysical quintessential type of matter field a priori to achieve the late 
acceleration. 

[Hi) It is observed that the instant of flip is delayed with dimensions, which is also 
supported by current observations. In this respect the HD accelerating models have an 
edge over the 4D ones. 

{iv) The most important thing in the second case, in our opinion, is the finding that 
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depending on some initial conditions, the effective equation of state during late evolution 
interpolates among ACDM, QCDM and Phantom type of expansion. In this respect our 
work recovers the effective equation of state (for large scale factor) for a recent work of 
Guo et al where a very generalised Chaplygin type of gas is taken. 

{v) Our solutions are general in nature because all the known results of 4D cosmology 
are recovered when d = 0. 
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